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Abstract: We experimentally investigate the non-Gaussian features of the 
phase-randomized coherent states, a class of states exploited in communica- 
tion channels and in decoy state-based quantum key distribution protocols. 
In particular, we reconstruct their phase-insensitive Wigner functions and 
quantify their non-Gaussianity. The measurements are performed in the 
mesoscopic photon-number domain by means of a direct detection scheme 
involving linear detectors. 
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1. Introduction 

Recently, weak phase-randomized or phase-averaged coherent states (PHAVs) have been suc- 
cessfully exploited to implement quantum key distribution (QKD) |T||2][3). It has been demon- 
strated that phase randomization is useful to improve the security of QKD, pointing out the 
relevance of a careful characterization of the employed sources 0E). With respect to coherent 
states, which are described by Gaussian Wigner functions (SJ, PHAVs clearly exhibit non- 
Gaussian features. Thus, the systematic study of the nature of PHAVs and the possibility to 
manipulate it can be considered of real interest in enhancing the performances of the commu- 
nication protocols in which such states are used (7). The complete characterization of optical 
states can be obtained by reconstructing the Wigner function. In practice, two alternative de- 
tection techniques, both based on interferometric schemes, can be used to achieve this goal. 
The typical method employed so far is the time-resolved optical homodyne tomography (see 
Ref. (HI for a review). Nevertheless, on the basis of the results derived in Ref. [9], it has been 
demonstrated that also direct detection can be used to obtain the Wigner function starting from 
photon-number distributions ifTQllTTIl . 

In this paper we report on the controlled generation of PHAVs Ifl2l and on their full char- 
acterization in terms of Wigner functions by exploiting a direct detection scheme involving 
linear detectors. In particular, we present a systematic study of their non-Gaussianity (nonG) 
in different experimental conditions and for different choices of the parameters that determine 
the states. Moreover, we demonstrate that it is possible to manipulate the nonG by linear pas- 
sive devices (beam splitters). In order to quantify the nonG of the states we adopt a recently 
proposed measure given in terms of the von Neumann entropy lfT3l[T4l . This investigation is 
performed in the mesoscopic photon-number regime as in this domain the optical states are 
robust with respect to losses and thus it is possible to better understand the influence of the 
different parameters on the nature of PHAVs. 

The paper is structured as follows. In Section [2] we introduce the concepts and formalism 
used throughout the paper. In particular we write the theoretical Wigner functions of PHAVs 
and of the state obtained by mixing two PHAVs at a beam splitter. Section [3] describes our 
experimental setup for the generation, manipulation and characterization of PHAVs and con- 
tains the experimental reconstruction of the Wigner functions of the generated states and the 
results we obtained for the nonG. We close the paper by drawing some concluding remarks in 
SectionlH 

2. Theory 

The density matrix pp of a single-mode PHAV of amplitude j3 is given by: 



where |j3) is a coherent state, j3 = |j3 1 e'^, and {|«)}, n = 0, 1,2, . . ., is the photon number basis. 
The corresponding Wigner function is 03]: 




(1) 




(2) 



Iq(z) being the modified Bessel function. The representation in terms of Wigner function em- 
phasizes both the non-Gaussian nature of PHAV and its phase-insensitive nature (i.e. rotational 
invariance about the origin of the phase space). 

In order to assess the nonG of the state pp, we consider the following measure fl3l: 

8[pp}=S(a)-S(pp), (3) 

where S(p) = — Tr[plnp] is the von Neumann entropy of the state p and a is a reference 
Gaussian state with the same covariance matrix as the state p . As the PHAV is a diagonal state 
in the photon number basis, its reference state is a thermal equilibrium state, with the same 
mean number of photons TV = |/3 1 2 f 13 1. In this case Eq. (01 reduces to: 

5 [p p ] = {N + 1 ) ln(N + 1 ) - NlnN + £ p„ In p n , (4) 

n 

where p„ = exp (-|j3 1 2 ) |/3 \ 2n /n\. 

Another interesting state we will investigate is that obtained from the interference of two 
PHAVs pp and pp at a beam splitter (BS) with transmissivity T, which can be useful for appli- 
cation to passive decoy state QKD [7]. In this case, if we consider only one of the two modes 
outgoing the beam splitter, we still obtain a state with a diagonal density matrix whose Wigner 
function can be written as follows J6) : 

W 2 _phav («; P\£\t)= / -^WpHAv(a-j3Vr^;/3v^) (5) 
Jo 2n 

where /3 = |j3|e"' > and the function in the integral is given by Eq. (0. Obviously, the Wigner 
function of the other outgoing mode can be obtained by replacing T with the reflectivity (1 — t). 
In the following we will refer to this state as 2-PHAV and, since it is diagonal, its nonG is still 
given by Eq. where now p n is the 2-PHAV non-trivial photon statistics reported in lfl2l . 



3. Experimental results 

The measurements we performed were aimed at reconstructing the Wigner functions of PHAVs 
and 2-PHAVs and at comparing their nonG for different choices of the experimental parame- 
ters. The generation of the class of the PHAVs was achieved by exploiting the second harmonics 
(@ 523 nm, 5-ps pulses) of a mode-locked Nd:YLF laser amplified at 500 Hz (High-Q Laser 
Production). According to the experimental setup sketched in Fig. [T] the single PHAV was 
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Fig. 1. Experimental setup. F: variable neutral density filter; BS: 50/50 beam splitter; Pz: 
piezoelectric movement; MF: multimode fiber (600 pm core). 



obtained by sending the light pulses to a mirror mounted on a piezo-electric movement. Its dis- 
placement, which was controlled by a function generator, was operated at a frequency of 100 Hz 
and covered a 12 jJ.m span. The 2-PHAV was produced by mixing at a BS two single PHAVs 
obtained by duplicating the phase-randomization system. A filter inserted in the path of one of 
the two PHAVs allowed us to change the balancing between the two fields. We implemented 
a direct detection scheme involving a photon-counting detector, namely a hybrid photodetec- 
tor (HPD, R10467U-40, maximum quantum efficiency ~0.5 at 500 nm, 1.4 ns response time, 
Hamamatsu), that is a detector not only endowed with a partial photon-counting capability, but 
also characterized by a linear response up to 100 photons. Thanks to its features, the HPD can 
actually operate in the mesoscopic domain, where the states are robust with respect to losses. 
The output of the detector was amplified (preamplifier A250 plus amplifier A275, Amptek), 
synchronously integrated (SGI, SR250, Stanford) and digitized (AT-MIO-16E-1, National In- 
struments). The estimated overall quantum efficiency of the detection chain is equal to 0.45. 
By using the strategy described in Refs. |[T6l[l7l . we can reconstruct the detected-photon dis- 
tributions that, as we have already proven in [|T5| , represent the basic elements to retrieve the 
Wigner function. 

The experimental reconstruction of this function was achieved by mixing at a BS the state 
under investigation, in the present case either the single PHAV or the 2-PHAV, with a coherent 
probe field whose amplitude and phase could be continuously varied. As both the states to be 
characterized are phase-insensitive, we reconstructed only a section of the Wigner distribution 
for fixed phase. In Fig. [2] we plot the experimental data (black dots + error bars) of a single 
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Fig. 2. Experimental reconstruction of a section of the phase-insensitive Wigner function 
of a PHAV (left panel), with |j8| 2 = 1.97, and of a 2-PHAV (right panel), with |j3| 2 = 1.03 
and |/3| 2 = 0.91. Black dots: experimental data; red line: theoretical expectations given in 
Eqs. ® and l[7}, respectively. 



PHAV (left panel) and of a balanced 2-PHAV (right panel), endowed with nearly the same 
mean number of detected photons, Mj = 1.97 and Mj = 1.94, respectively. In each panel we 
also show the theoretical expectations (red line) for the PHAV and 2-PHAV, respectively 031 : 

WpHAv(v / ?«)=WpHAv(v / ?«K^ ?(|a| + l/3|) , (6) 

^ being the overall (spatial and temporal) overlap between the probe and the PHAV, and: 

W 2 _phav(v / ¥«) = W2-PHAv(v / ?p«)^ [ ^ Tr¥|a| + ^ /Tr¥(l,3| + l/3|)1 , (7) 

where £,p describes the overall overlap between the probe and the 2-PHAV and £,s the overall 
overlap between the two components of the 2-PHAV. In both the Eqs. ^ and (0, |/3 1 2 and |j3 1 2 



are now the mean numbers of photons we measured (see Fig. [2]), thus including the quantum 
efficiency. In fact, it is worth noting that for classical states the functional form of the Wigner 
function is preserved also in the presence of losses and its expression, given in terms of detected 
photons, reads W(a) = ^ Ym=o(~ l) m Pm'a' where Pma represent the detected-photon distribu- 
tions of the state to be measured displaced by the probe field |[T5l . As testified by the very high 
values of the overlaps £, , £,p an d 4s reported in the panels of Fig. [2 we actually achieved a very 
good superposition in aligning the system. In our opinion, having overcome this criticality can 
be considered a demonstration of the robustness of our setup and also a fundamental step on 
the way to the exploitation of these states for possible applications. 

We characterized the nonG amount of the reconstructed states by considering the measure 
in Eq. @, which, however, is actually given in terms of photons. As we have not directly 
access to the number of photons, we calculated a similar expression, e, for detected pho- 
tons, which represents a lower bound to nonG USUI)- We found £phav = 0.207 ± 0.004 and 
£2-phav = 0.036 ± 0.005. From this comparison it emerges that a PHAV is more non-Gaussian 
than the 2-PHAV, thus demonstrating that a Wigner function exhibiting a dip in the origin of 
the phase space is more non-Gaussian than one characterized by a peak in the origin followed 
by a "shoulder". The results prove that combining two non-Gaussian states does not neces- 
sarily lead to an increase of the overall nonG. Indeed, to better understand the nature of the 
"shoulder", we generated balanced 2-PHAVs at different mean values and reconstructed their 
Wigner function. The experimental data are shown in Fig.[3]together with the theoretical curve. 
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Fig. 3. Experimental reconstruction of a section of the phase-insensitive Wigner function of 
a 2-PHAV at different mean values Mj (see the text for details). Black dots: experimental 
data; red line: theoretical expectations given in Eq. 10. 



We can observe that the "shoulder" is more and more evident at increasing the mean number 
of photons of the balanced 2-PHAV. Indeed, the calculation of the nonG amount is in accor- 
dance with this result, as it is an increasing function of the mean number of photons (left panel 
£2-phav = 0.03 1 ± 0.005, Mj = 1.61, and right panel £ 2 -phav = 0.061 ± 0.006, M T = 2.92). 
For the sake of clarity, in the left panel of Fig. [4] we show the behavior of £ as a function of 
the mean total energy of the 2-PHAV. To achieve a more complete characterization of this class 
of states, we also studied the behavior of the nonG amount by varying the balancing between 
the two components of the 2-PHAV and keeping the mean number of detected photons of the 
overall state fixed. The lower bound £ of the nonG amount as a function of the ratio between 
the two components is plotted in the right panel of Fig.|4] as one may expect, it monotonically 
decreases at increasing the balancing. In fact, the most unbalanced condition reduces to the 
case in which there is only a single PHAV, whereas the most balanced one corresponds to have 
a balanced 2-PHAV. In Fig. |4]we note that the experimental values of £ are larger than those 
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Fig. 4. Left panel: nonG measure e for the detected number of photons as a function of 
the mean number of photons of balanced 2-PHAVs (see the text for details). Right panel: £ 
as a function of the balancing between the two components of the 2-PHAV, at fixed mean 
number of detected photons Mj = Mi +M2 of the overall state. Black dots: experimental 
data; red squares: theoretical expectation. 



expected from theory and the difference increases at high mean numbers (left panel) or at high 
unbalancing (right panel). We ascribe this behavior to small errors (few percent) in the determi- 
nation of the mean values of detected photons that can occur at relatively high intensities due 
to some saturation of the detector. Indeed, the sensitivity of £ to the mean number of photons 
is testified by the error bars associated to the experimental data. This critical character suggests 
the exploitation of such nonG measure to test decoy states based on 2-PHAVs Q. 

4. Conclusion 

In conclusion, we have presented an experimental investigation of the non-Gaussian nature of 
the class of PHAVs by reconstructing the Wigner function of the states in several different 
conditions. We have also quantified the nonG amount by exploiting a nonG measure defined in 
terms of von Neumann entropy and proved the consistency of the two approaches. It is worth 
noting that the characterization and manipulation of the nonG could be a useful tool for the 
enhancement of the QKD protocols in which PHAVs are employed. 

Moreover, the very good agreement between experimental data and theory reinforces the 
possibility of using this class of states for further investigations on the nonG of optical states 
and in particular for the realization of more complex experimental schemes, also involving 
nonclassicality. Indeed, the good quality of the results makes our system particularly useful for 
quantum state engineering and quantum information processing, despite the many criticalities 
and difficulties in obtaining a good randomization of the phases, in improving the superposition 
of more than two light fields (such as in the case of the Wigner function reconstruction of the 
2-PHAV) and in maintaining the alignment for hours and for different kinds of measurements. 

Finally, it is worth noting that the measurements were actually performed in the mesoscopic 
photon-number domain, as the reconstruction of the Wigner functions of the states under in- 
vestigation was achieved by displacing them with a coherent field whose intensity was changed 
from zero up to four times the mean value of the states themselves. 
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